Angle-Dependent Spin- Wave Resonance Spectroscopy of (Ga,Mn)As Films 
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A modeling approach for standing spin-wave resonances based on a finite-difference formulation 
of the Landau-Lifshitz-Gilbert equation is presented. In contrast to a previous study [Bihler et ai, 
Phys. Rev. B 79, 045205 (2009)], this formalism accounts for elliptical magnetization precession and 
magnetic properties arbitrarily varying across the layer thickness, including the magnetic anisotropy 
parameters, the exchange stiffness, the Gilbert damping, and the saturation magnetization. To 
demonstrate the usefulness of our modeling approach, we experimentally study a set of (Ga,Mn)As 
samples grown by low-temperature molecular-beam epitaxy by means of angle-dependent stand- 
ing spin-wave resonance spectroscopy and electrochemical capacitance-voltage measurements. By 
applying our modeling approach, the angle dependence of the spin-wave resonance data can be re- 
produced in a simulation with one set of simulation parameters for all external field orientations. 
We find that the approximately linear gradient in the out-of-plane magnetic anisotropy is related 
to a linear gradient in the hole concentrations of the samples. 

PACS numbers: 75.50.Pp, 76.50.+g, 75.70.-i, 75.30.Ds 
Keywords: (Ga,Mn)As; spin wave resonance; magnetic anisotropy 
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I. INTRODUCTION 

Due to their particular magnetic properties, in- 
cluding magnetic anisotropy^— anisotropic magneto- 
resistance^ and magneto-thcrmopowerj 6 - in the past 
years ferromagnetic semiconductors have continued to 
be of great scientific interest in exploring new physics 
and conceptual spintronic devices The most promi- 
nent ferromagnetic semiconductor is (Ga,Mn) As, where a 
small percentage of Mn atoms on Ga sites introduces lo- 
calized magnetic moments as well as itinerant holes which 
mediate the ferromagnetic interaction of the Mn spins (p- 
d exchange interaction)^ Both theoretical and experi- 
mental studies have shown that the magnetic anisotropy, 
i.e., the dependence of the free energy of the ferromagnet 
on the magnetization orientation, depends on the elas- 
tic strain and the hole concentration in the (Ga,Mn)As 
layer ; 12 ' 13 opening up several pathways to manipulate the 
magnetic anisotropy of (Ga,Mn)AsJ^~— 

A common spectroscopic method to probe the mag- 
netic anisotropy of ferromagnets and in particular 
(Ga,Mn)As, is angle-dependent ferromagnetic resonance 
(FMR)^ - — where FMR spectra are taken as a func- 
tion of the orientation of the external magnetic field. If 
the magnetic properties of the ferromagnet are homo- 
geneous, a zero wave vector (k = 0) mode of collec- 
tively, uniformally precessing magnetic moments couples 
to the microwave magnetic field, e.g., in a microwave cav- 
ity, allowing for a detection of the magnetization preces- 
sion. The resonance field of this mode, referred to as 
uniform resonance magnetic field, depends on the em- 
ployed microwave frequency and the magnetic anisotropy 



parameters. Thus, by recording FMR spectra at differ- 
ent orientations of the external field with respect to the 
crystal axes, the anisotropy parameters can be deduced 
from the experiment. However, if the magnetic prop- 
erties of a ferromagnetic layer are non-homogeneous or 
the spins at the surface and interface of the layer are 
pinned, non-propagating modes with k =/= 0, referred to 
as standing spin- wave resonances (SWR), can be excited 
by the cavity field and thus be detected in an FMR ex- 
periment. On one hand this can hamper the derivation 
of anisotropy parameters, on the other hand a detailed 
analysis of these modes can elucidate the anisotropy pro- 
file of the layer and the nature of spin pinning condi- 
tions. Furthermore, the excitation of spin waves is of 
topical interest in combination with spin-pumping j 2 ^— 
i.e., the generation of pure spin currents by a precessing 
magnetizationi^r— In this context, the exact knowledge 
of the magnetization precession amplitude as a function 
of the position coordinate within the ferromagnet is of 
particular importance^ 

Several publications report on SWR modes in 
(Ga,Mn)As with a mode spacing deviating from what is 
expected according to the Kittel model for magnetically 
homogeneous films with pinned spins at the surface^ - — 
These results have been attributed to an out-of-plane 
anisotropy field linearl y 31 : 36 or quadratically varying^ 3 . - — 
as a function of the depth into the layer, as well as to 
specific spin pinning conditions at the surface and at the 
interface to the substrate^ While most of these studies 
have focused on the spacings of the resonance fields when 
modeling SWR measurements, in Ref. [36| a more sophis- 
ticated approach, based on a normal mode analysis ) 37 : 38 
was employed to model resonance fields as well as rela- 



tive mode intensities for the external field oriented along 
high-symmetry directions, assuming a circularly process- 
ing magnetization. 

In this work, we present a more general modeling ap- 
proach for SWR, based on a finite-difference formulation 
of the Landau-Lifshitz-Gilbert (LLG) equation. This ap- 
proach holds for any orientation of the external mag- 
netic field and accounts for elliptical magnetization pre- 
cession [Sec. QTJ. It allows for a simulation of arbitrar- 
ily varying profiles of the magnetic properties across the 
thickness of the film, including vatiations of the mag- 
netic anisotropy parameters, the exchange stiffness, and 
the Gilbert damping parameter. As the result of the sim- 
ulation, we obtain the Polder susceptibility tensor as a 
function of the depth within the ferromagnet. Based on 
this result, the absorbed power upon spin wave resonance 
and the magnetization precession amplitude as a function 
of the depth can be calculated for any orientation of the 
external magnetic field. 

We apply our modeling approach to a set of four 
(Ga,Mn)As samples epitaxially grown with different 
V/III flux ratios [Sec. IIII] . motivated by the obser- 
vation that V/III flux ratios of < 3 lead to a gra- 
dient in the hole concentration p [Ref. HH, which in 
turn is expected to cause non-homogeneous magnetic 
anisotropy parameters . 31 ! 36 Electrochemical capacitance- 
voltage (ECV) measurements revealed a nearly linear 
gradient in p across the thickness of the layers investi- 
gated. To show that our modeling approach is capa- 
ble of simulating SWR spectra for arbitrary magnetic 
field orientations, angle-dependent SWR data were taken 
and compared with the model using one set of magnetic 
parameters for each sample, revealing gradients in the 
uniform resonance magnetic fields. We discuss the in- 
fluence of the gradient in p on the observed uniform 
resonance field gradients as well as possible influences 
of strain and saturation magnetization gradients on the 
observed out-of-plane anisotropy profile. It should be 
emphasized, however, that the objective of this work is 
to show the usefulness of our modeling approach, while 
a detailed investigation of the origin of the gradient in 
the out-of-plane magnetic anisotropy profile and there- 
fore a detailed understanding of the particular materials 
physics of (Ga,Mn)As is beyond the scope of this study. 
Finally, we summarize our results and discuss further po- 
tential applications of this work [Sec. HVj . 



II. THEORETICAL CONSIDERATIONS 

In this section, we provide the theoretical framework 
necessary to describe the full angle dependence of the 
spin- wave resonance spectra presented in Sec. IIIII Refer- 
ring to the coordinate system depicted in Fig. [TJ we start 
from the canonical expression for the free enthalpy den- 
sity (normalized to the saturation magnetization M) for 
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FIG. 1: (color online) Relation between the two coordinate 
systems employed. The (a;, y, z) frame of reference is spanned 
by the cubic crystal axes, while the (1,2,3) coordinate sys- 
tem is determined by the equilibrium orientation of the mag- 
netization (3-direction) and two transverse directions, the 2- 
direction being parallel to the film plane; the latter system is 
z and fioH dependent, as described in the text. 

a tetragonally distorted (Ga,Mn)As filmi^i&ii 
G = const — hqH ■ m + B aa im 2 z + B±\_m z 

+ B 4\\( m t +™$) + ^-BlloOa: - m v) 2 - ( l ) 

Here, [AqH is a static external magnetic field, Booi 
is a uniaxial out-of-plane anisotropy parameter, re- 
flecting shape and second-order crystalline anisotropy, 13 
£?4j_, B411, and B x i are fourth-order crystalline and 
second-order uniaxial in-plane anisotropy parameters, 
respectively m x , m y , m z denote the components of the 
normalized magnetization vector m(z) = M(z)/M(z) 
along the cubic axes [100], [010], and [001], respectively. 
We assume the magnetic properties of the layer to be ho- 
mogeneous laterally (within the xy plane) and inhomoge- 
neous vertically (along z); the anisotropy parameters in 
Eq. {TJ and the magnetization are consequently a func- 
tion of the spatial variable z. To obtain the anisotropy 
parameters from Eq. (TTJ in units of energy density, it 
would therefore be required to know the z dependence 
and the absolute value of M. 

The minimum of Eq. (JTJ) determines the equilibrium 
orientation of the magnetization, given by the angles 
6*o = (z) and 0o = 4>o( z ): cf. Fig. [TJ To describe the 
magnetization dynamics, we introduce a new frame of 
reference (1,2,3) shown in Fig. [TJ in which the equilib- 
rium orientation of the magnetization mo coincides with 
the axis 3. For small perturbations, the magnetization 
precesses around its equilibrium with finite transverse 
components of the magnetization rrii (i — 1,2) as illus- 
trated in the inset in Fig. [TJ The transformation between 
the two coordinate systems is given in the Appendix [A] 
by Eqs. (|A1[) and (|A2[) . We write for the (normalized) 
magnetization 



3 




0{m\, m 2 ) 



(2) 



The evolution of the magnetization under the influence 
of an effective magnetic field ^o-Hcff is described by the 
LLG equatior^ 3 - 



(3) 



d t m = —7m x fj, H c ff + am x d t m, 



where 7 is the gyromagnetic ratio and a a phcnomcno- 
logical damping parameter. The effective magnetic field 
is given by 3 - 6 - 



fi H eS = -V m G + -r|V 2 M + fi h{t) 



M 



(4) 



where V m = (d m i, 9 m 2, d m3 ) is the vector differential op- 
erator with respect to the components of ra, D s = 2A/M 
is the exchange stiffness with the exchange constant A, 
V 2 is the spatial differential operator V 2 = 9 2 + <9 2 + <9 2 , 
and h(t) — hoe~ luJt is the externally applied microwave 
magnetic field with the angular frequency oj; h(t) is ori- 
ented perpendicularly to (IqH. Since the magnetic prop- 
erties are independent of x and y, Eq. ([3]) simplifies to 



d t m = —7m x [-V m G + D s m" + fi h(t)] 



■ am x dtm, 
(5) 

with m" — <9 2 m, neglecting terms of the order of m 2 (for 
i = 1, 2). By definition of the (1, 2, 3) coordinate system, 
the only non-vanishing component of V m G in the equi- 
librium is along the 3-direction. For small deviations of 
m from the equilibrium we find 44 



V m G = 



Gum,! + G 2 \m 2 
G12TO1 + G 22 m 2 
G 3 



(6) 



abbreviations Gi = 
<i=m ; the explicit ex- 



where we have introduced the 
d mi G\ m—rno 

and Gij = d m% d mj G\ 
pressions for these derivatives are given in the Appendix 

m 

In the following, we calculate the transverse magnetiza- 
tion components assuming a harmonic time dependence 



The linearized LLG equation, consider- 



ing only the transverse components, reads as 



Hu H 12 \ / mi 

H 2 1 H 2 2 J V m 2 



= Mo 



hi 



(7) 



where we have introduced the abbreviations Hu = Gu — 
G3 - iaw/j, ilia = H|i = G12 + «w/7, and H 2 2 = G22 - 
G3 — iacj/7. We have dropped all terms which are non- 
linear in mi and products of m.j with the driving field. 

Resonant uniform precession of the magnetization 
(to" = 0) occurs at the so called uniform resonance field 



Mo-ffuni(z), which is found by solving the homogeneous 
(h = 0) equation 

H 11 (z)H 22 (z) - H 12 (z)H 2 i(z) = 

^ (G n - G 3 )(G 2 2 - G 3 ) - G 2 12 = (~) (8) 

for hqH, neglecting the Gilbert damping (a = 0). Equa- 
tion (JSJ can be used to derive anisotropy parameters from 
angle-dependent FMR spectra. As extensively discussed 
by Baselgia et al.~, using Eq. (jHJ is equivalent to using 
the method of Smit and Beljers, which employs second 
derivatives of the free enthalpy with respect to the spher- 
ical coordinates! 41 ! 45 ' 46 

To illustrate the role of the uniform resonance field in 
the context of spin- wave resonances, we consider the spe- 
cial case where magnetization is aligned along the [001] 
crystal axis (8q — 0), before we deal with the general 
case of arbitrary field orientations. Neglecting the uniax- 
ial in-plane anisotropy (-B110 = 0) since this anisotropy is 
typically weaker than all other anisotropics ; 13 ' 41 we find 
G 3 = -fi H + 2B m i + AB i± _ and Gu = G 22 = G 12 = 0, 
resulting in the uniform resonance field 



/xoiCr («) - w/7 + 2B 001 (z) + 4S 4± (z). 



(9) 



To find the eigenmodes of the system, we consider the 
unperturbed and undamped case, i.e., a — and h = 
in Eq. ([7]). With to-2 = imi = fh we find the spin- wave 
equation 



D s m" + ^ H l Zi(z)m - VoHm 



(10) 



in agreement with Ref . l36l . The relation of the anisotropy 
parameters defined in Ref. to the ones used here 
is given by B 001 = K™°/M + B ll0 , B ll0 = -K^/M, 

2B 4± = -K^-JM, and 2B 4 || = -K^/M. Equation (JTOj) 
is mathematically equivalent to the one-dimensional 
time-independent Schrodinger equation, where the 
uniform resonance field corresponds to the potential, 
to to the wave function, [j,qH to the energy, and D s is 
proportional to the inverse mass. To calculate the actual 
precession amplitude of the magnetization, the coupling 
of the eigenmodes of Eq. (| 10[) to the driving field is 
relevant, which is proportional to the net magnetic 
moment of the mode . 36 ' 38 In analogy to a particle in a 
box, the geometry of the uniform resonance field as well 
as the boundary conditions determine the resonance 
fields and the spatial form of the precession amplitude. 
For the remainder of this work, we assume the spins 
to exhibit natural freedom at the boundaries of the 
film, i.e., d z rh — fh' — Q at the interfaces ) 36 ' 47 since 
these boundary conditions have been shown to describe 
the out-of-plane SWR data of similar samples well.— 
To graphically illustrate the influence of the uniform 
resonance field on the SWR modes, we consider in Fig. [2] 
a ferromagnetic layer with a thickness of 50 nm with 
constant magnetic properties across the layer (a) and 
with a linearly varying uniform resonance field (c); in 
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FIG. 2: Simulation to demonstrate the influence of the uni- 
form resonance field fioH^J on the SWR modes for mo 1 1 [001] , 
assuming circular precession. In (a), fto-ffuni i s se ^ to be con- 
stant across the layer, while in (c) it varies linearly (blue, 
dashed lines), in analogy to a square potential and a trian- 
gular potential, respectively. The dotted black lines are the 
resonance fields, calculated assuming boundary conditions of 
natural freedom, see text. The solid red lines show the eigen- 
modes of the system, i.e., the precession amplitude m of the 
magnetization; for each mode the dotted line corresponds to 
fh = 0. As can be seen in (a), for a constant uniform reso- 
nance field the first mode occurs at the uniform resonance field 
and exhibits a constant precession amplitude across the layer, 
i.e., an FMR mode. The second and third mode (higher-order 
modes are not shown) exhibit a non-uniform magnetization 
profile. In order to couple to the driving field the modes need 
to have a finite net magnetic moment. As can be seen in 
(a), the positive and negative areas of the second and third 
mode are equal, thus these modes are not visible in the SWR 
spectrum (b). This is in contrast to the case of the linearly 
varying uniform resonance field (c) where the mode profile is 
given by Airy functions, which have a nonzero net magnetic 
moment also for the second and third mode, resulting in a 
finite SWR intensity of these modes (d). The spectra in (b) 
and (d) were calculated by integrating over the eigenmodes 
m and convoluting the square of the result with Lorentzians. 



both cases we assume D s — 13 Tnm~ 2 , a similar value 
as obtained in previous studies^ For these conditions, 
we numerically solve Eq. (|10[) by the finite difference 
method described in the Appendix lB 11 in order to obtain 
the resonance fields (eigenvalues) and the z dependence 
of the transverse magnetic moments (eigenfunctions) . 
To which amount a mode couples to the driving field is 
determined by the net magnetic moment of the mode, 
which is found by integrating rh{z) over the thickness of 
the film. For the magnetically homogeneous layer, the 
only mode that couples to the driving field is the uniform 



precession mode at noH^j, since modes of higher order 
have a zero net magnetic moment [Fig. [2] (a)], resulting 
in one resonance at the uniform resonance field, cf. Fig. [2] 
(b). For the non- uniform layer, with [IqH®®1(z) linearly 
varying across the film, the mode profile is given by Airy 
function o 31 ^ 36 ' 38 and various non-uniform modes couple 
to the driving field, resulting in several spin-wave reso- 
nances with their amplitude proportional to the square 
of the net magnetic momen t 36 ' 38 of the corresponding 
mode, cf. Fig. H(c) and (d). 

We now turn to the general case of arbitrary field ori- 
entations. Due to the magnetic anisotropy profile, the 
magnetization orientation is a priori unknown and a func- 
tion of z and [j.qH. Furthermore, the assumption of a 
circularly precessing magnetization is not generally jus- 
tified. To solve Eq. (J7J for arbitrary field orientations, 
we employ a finite difference method as outlined in the 
Appendix lB 21 By solving Eq. ([7]), we obtain the z depen- 
dent generalized Polder susceptibility tensor x{lM)H,z), 
which relates the transverse magnetization components 
Mj(z) = M(z)m,i(z) with the components of the driving 
field by 



M 1 
M 2 



hi 
h 2 



(11) 



In a microwave absorption measurement, the components 
Mi which are out-of-phase with the driving field are de- 
tected. The absorbed power density is related to the 
imaginary part of x([XoH, z) and can be calculated by 48 



P = 



2z Q 



Im 



K,h-2 ) x(jM)H,z) 



hi 
h 2 



(12) 

where zq is the thickness of the ferromagnetic layer. Note 
that the position coordinate z is negative in the film, 
cf. Fig. ffl 

To obtain an impression of how gradients in differ- 
ent anisotropy parameters influence the SWR spectra, 
we plot in Fig. [3] simulated SWR spectra together with 
the magnetization precession cone as a function of depth 
in the ferromagnetic layer. We assume a constant sat- 
uration magnetization (its value is not relevant for the 
outcome of the simulation), a constant exchange stiff- 
ness D s — 35 Tnm 2 unless otherwise specified, a = 0.09, 
and i?ooi = 90 mT, £> 4 || = —50 mT, = 15 mT. 

In Fig. [3] (a), we assume -Boot to vary across the layer 
thickness according to Bqqi(z) = £?ooi — &ooi x z with 
&ooi = —0.8 mT/nm. Figure [3] (a i) shows the simu- 
lated SWR spectra calculated by taking the first deriva- 
tive of Eq. (|T2|) with respect to /j,qH for different angles 
tfj defined in the inset in Fig. [3] (c iv). We observe sev- 
eral SWR modes for /io-ff||[001] which become less as 
/j,qH is tilted away from [001]. At tp = 40° only one 
mode is visible while for ip — 0° we again observe mul- 
tiple SWR modes. This observation can be understood 
by considering the uniform resonance fields as a func- 
tion of the depth for these orientations. In Fig. [3] (a 
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FIG. 3: Atlas illustrating the influence of gradients in the anisotropy parameters on SWR spectra. In (a) all anisotropy 
parameters are kept constant with the values given in the text, except Booi which is varied linearly. Correspondingly, in (b) 
and (c) B4± and Bm\ were varied linearly, respectively. Panels (i) show the first derivative of simulations using Eq. (|12p with 
respect to fioH and panels (ii)-(iv) show the precession cone Im(m*m2 — mimjj) in a color plot together with the uniform 
resonance field /j.oH un i(z) (dashed blue lines) at three different external field orientations; the black dotted lines indicate the 
resonance field positions of the modes. Panel (a i) additionally shows the influence of a linear gradient in the exchange stiffness 
parameter on the spin-wave spectra, see text for further details and discussion. 



ii) -(a iv), we show the uniform resonance field (dashed 
blue line) for ip — 0°, ip — 30°, and ifi = 90°, respec- 
tively, together with the magnetization precession cone 
Im(m 1 m2 — mim^) in a contour plot as a function of 
depth and fioH. At ip = 90°, the uniform resonance field 
varies strongly across the film, which can be understood 
by considering Eq. ([§]). This results in several spin wave 
modes with their resonance fields indicated by dotted 
lines. 

For other field orientations, the formula for the uni- 
form resonance field can also be derived but results in a 
longer, more complex equation than Eq. @. Important 
in this context is that positive values of -Boor lead to an 
increase (decrease) of the resonance field for the magneti- 
zation oriented perpendicular (parallel) to the film plane, 
accounting for the reversed sign of the slopes of ^aH unl in 
Fig. [3] (a ii) and (a iv). Consequently, in between those 
two extreme cases fioH un i must be constant across the 
layer for some field orientation, in our case for tp = 30°, 
resulting in a single SWR mode, cf. Fig. [3] (a i) and (a 

iii) . In addition to the SWR simulations with constant 



D s , we plot in Fig. [3] (a i) simulated SWR spectra with D s 
varying linearly across the film with D s = 35 — 65 Tnm 2 
(blue, dotted lines) and D s = 35 — 5 Tnm 2 (green, dotted 
lines). A decreasing D s leads to a decreasing spacing in 
the modes and vice versa for an increasing D s as can be 
seen, e.g., for // J?||[001]. 

In Fig. [3] (b), we consider the case where all magnetic 
parameters are constant with the values given above, ex- 
cept -B4_l(z) = i?4_L — 64^ x z with 64^ = —0.4 mT/nm. 
As evident from Eq. ([9]), this results in the same slope of 
Mo-ffuni for i\) = 90° as in the case above where we varied 
£?ooi only, cf. Fig. [3J (a iv) and (b iv). In contrast to the 
case depicted in (a), however, here for -0 = 0° the uniform 
resonance field is constant. This can be understood when 
evaluating the parameters that enter in the calculation of 
the uniform resonance field [Eq. ©]. If m is in the film 
plane, none of the parameters in Eqs. (|A4[) - (|A6[) depends 
on i?4j_, resulting in a constant uniform resonance field 
for -0 = 0°. As m is tilted away from the film plane, 
i?4_i_ enters in some of the terms Eqs. (|A4|) - (| A6|) . As a 
consequence, /io-ff U ni varies, first such that it increases 
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[cf. Fig. [3] (b iii)] and finally, such that it decreases as a 
function of depth [cf. Fig. (b iv)]. 

Finally, we discuss the case where all parameters are 
constant except B4\\(z) = B±\\ — 6411 x z with 6411 = 
—0.4 mT/nm [Fig. |5] (c)]. Here, fj,oH unl is constant for 
ip = 90° as predicted by Eq. © . As m is tilted away from 
[001] a varying Bm\ leads to a varying uniform resonance 
field as shown in Fig. [3] (c ii) and (c iii). Here, a sign 
reversal of the slope as it was the case in Fig. [3] (a) and 
(b) does not take place and multiple resonances occur, 
starting from ip = 60° [Fig. [3] (c i)] . 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

(Ga,Mn)As samples with a nominal Mn concentration 
of w 4% were grown on (OOl)-oriented GaAs substrates 
by low-temperature molecular-beam epitaxy at a sub- 
strate temperature of 220°C using V/III flux ratios of 
1.1, 1.3, 1.5, and 3.5, referred to as sample A, B, C, and 
D, respectively. The layer thickness was 210-280 nm as 
determined from the ECV measurements, cf. Fig. [4] For 
samples with V/III flux ratios <3a gradient in the hole 
concentration has been reported^ hence this set of sam- 
ples was chosen to study the influence of a gradient in p 
on the out-of-plane magnetic anisotropy. Further details 
on the sample growth can be found in Refs. [1^ and lU 

The hole concentration profile of the as-grown 
(Ga,Mn)As layers were determined by ECV profiling us- 
ing a BioRad PN4400 profiler with a 250 ml aqueous 
solution of 2.0 g NaOH+9.3 g EDTA as the electrolyte. 
For further details on the ECV analysis see Ref. H^. The 
results of the ECV measurements for the layers investi- 
gated are shown in Fig. 2] (a). Except for the sample with 
V/III=3.5, they reveal a nearly linearly varying hole con- 
centration across the layer thickness with different slopes 
and with the absolute value of the hole concentration 
at the surface of the layer varying by about 20%. The 
profiles are reproducible within an uncertainty of about 
15%. 

To investigate the magnetic anisotropy profiles of the 
samples, we performed cavity-based FMR measurements, 
using a Bruker ESP300 spectrometer operating at a mi- 
crowave frequency of 9.265 GHz (X-band) with a mi- 
crowave power of 2 mW at T = 5 K; we used magnetic 
field modulation at a frequency of 100 kHz and an ampli- 
tude of 3.2 mT. Since we are mainly interested in the out- 
of-plane magnetic anisotropy, we recorded spectra for ex- 
ternal magnetic field orientations within the crystal plane 
spanned by the [110] and [001] crystal axes in 5° steps, 
cf. the inset in Fig. [SJ For each orientation, the field was 
ramped to 1 T in order to saturate the magnetization 
and then swept from 650 mT to 250 mT; the spectra for 
the samples investigated are shown in Fig. [S] 

We start by discussing qualitative differences in the 
spectra. The samples A and B exhibit several pro- 
nounced resonances for the external field oriented along 




depth (nm) 

FIG. 4: (a) The hole concentration in the different (Ga,Mn)As 
samples is shown as a function of the depth within the layers 
as determined by ECV profiling, (b) The uniform resonance 
fields (ioH®nl(z) for the four samples obtained from the sim- 
ulations for the out-of-plane orientation of the external field 
(ip = 90°) as a function of the depth. 

[001], which we attribute to standing spin- wave reso- 
nances [Fig. [5] (a) and (b)]. For these samples, the [001] 
direction is the magnetically hardest axis since at this 
orientation the resonance field of the fundamental spin- 
wave mode is larger than at all other orientations. As the 
external field is rotated into the film plane, the resonance 
position of this mode gradually shifts to lower field val- 
ues as expected for a pronounced out-of-plane hard axis. 
In contrast, the samples C and D exhibit the largest res- 
onance fields for a field orientation of 50-60° [Fig. [5] (c) 
and (d)] pointing to an interplay of second- and fourth- 
order out-of-plane anisotropy with different signs of the 
corresponding anisotropy parameters. These samples ex- 
hibit spin-wave resonances as well, however they are less 
pronounced than for samples A and B. 

To quantitatively model the spin-wave spectra we nu- 
merically solve for each magnetic field orientation the 
spin- wave equation ([7]) by the finite difference method as 
outlined in the Appendix IB 21 Although this method al- 
lows for the modeling of the SWR for arbitrary profiles 
of the anisotropy parameters, the exchange stiffness, the 
Gilbert damping parameter, and the saturation magne- 
tization, we assume the parameters to vary linearly as 
a function of z. This approach is motivated by the lin- 
ear gradient in the hole concentration, which in first ap- 
proximation is assumed to cause a linear gradient in the 
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FIG. 5: The spin-wave resonance data (dotted, blue lines) are shown together with simulations (red, solid lines) using the 
numerical procedure described in the text and in the Appendix IB 21 The data were obtained as a function of the external 
magnetic field orientation and magnitude for samples with a V/III flux ratio of (a) 1.1, (b) 1.3, (c) 1.5, and (d) 3.5. The 
rotation angle ip is defined in the inset and the parameters used for the simulations are summarized in Tab. [I] 
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anisotropy parameters, resulting in the spin-wave reso- 
nances observed in the samplesj 31 ' 36 In Tab. HI we have 
summarized the parameters used in the simulation for 
the different samples. The parameters in capital let- 
ters denote the value at the surface of the sample while 
the ones in lower-case letters denote the slope of this 
parameter; e.g., the z dependence of the second-order, 
uniaxial out-of-plane anisotropy parameter is given by 
-Booi(z) = ^001 — frooi x z - The layer thickness used for 
the simulation can be inferred from Fig. 2] (a) and was de- 
termined from the ECV data under the assumption that 
at the position where the hole concentration rapidly de- 
creases the magnetic properties of the layer abruptly un- 
dergo a transition from ferromagnetic to paramagnetic. 
For the simulations, we divided each film into n = 100 
layers with constant magnetic properties within each 
layer. For the gyromagnetic ratio we used 7 = g^/h 
with g = 221 

As result of the simulation we obtain the Polder sus- 
ceptibility tensor xi^oH, z ) and the transverse magneti- 
zation components as a function of z and (j,qH. Addition- 
ally, we obtain the z dependence of the uniform resonance 
field by solving Eq. ((8]) for each field orientation. In an 
SWR absorption experiment with magnetic field mod- 
ulation, the obtained signal is proportional to the first 
derivative of the absorbed microwave power with respect 
to the magnetic field. Thus, we calculate the absorbed 
power using the simulated susceptibility and Eq. (|12[) and 
numerically differentiate the result in order to compare 
the simulated SWR spectra with the experiment. Addi- 
tionally, we use a global scaling factor, accounting, e.g., 
for the modulation amplitude, which is the same for all 
field orientations, and we multiply all the simulated data 
with this factor. In Fig.[5j we plot the experimental data 
together with the simulations using the parameters given 
in Tab. HI demonstrating that a reasonable agreement be- 
tween theory and experiment can be found with one set of 
simulation parameters for all magnetic field orientations 
for each sample. 

We will now exemplarily discuss the angle dependence 
of the SWR spectrum of sample A shown in Fig. [SJa) 
based on the uniform resonance field and the resulting 
magnetization mode profile obtained from the simulation. 
To this end, we plot in Fig. [5] (a)-(c) the magnetization 
precession amplitude Im(m|m2 — rri\m\) for selected ex- 
ternal field orientations as a function of depth and ex- 
ternal magnetic field in a contour plot, together with the 
corresponding uniform resonance field. In Fig. [5] (d)-(f), 
we show for each external field orientation a magnifica- 
tion of the corresponding SWR spectrum together with 
the simulation. Note that in contrast to the normal-mode 
approach (Appendix IB 1[) used to calculate the modes in 
Fig. [5J where the coupling of each mode to the cavity 
field has to be found by integration, the approach elabo- 
rated in the Appendix IB 21 directly yields the transverse 
magnetization components, already accounting for the 
coupling efficiency and the linewidth. Further, the ap- 
proach presented in the Appendix IB 21 is also valid when 



the difference in the resonance fields of two modes is com- 
parable with or smaller than their linewidth, in contrast 
to the normal-mode approach 3 ^. 

If the external field is parallel to the surface normal 
('tjj = 90°) the uniform resonance field varies by about 
350 mT across the film thickness [cf. the dashed line 
in Fig. [5] (a)], resulting in several well-resolved stand- 
ing spin-wave modes. The spin-wave resonance fields are 
plotted as dotted lines in Fig. [5] (a); since the spacing of 
the resonance fields is larger than the SWR linewidth, 
the modes are clearly resolved, cf. Fig. [5] (a) and (d). 
In the simulation two regions with different 6001 values 
were used in order to reproduce the spacing of the higher- 
order spin- wave modes found in the experiment. Using 
the same slope as in the first 100 nm for the entire layer 
would lead to a smaller spacing between the third and 
higher order modes. Instead of defining two regions with 
different slopes &ooij & gradient in the exchange stiffness 
with positive slope could also be used to model the exper- 
imentally found mode spacing as discussed in the context 
of Fig.[3J Since the exchange interaction in (Ga,Mn)As is 
mediated by holes^ and p decreases across the layer, we 
refrain from modeling our results with a positive gradient 
in D s . Further, the results in Ref. [36| rather point to a 
negative gradient in D s in a similar sample. However, a 
decreasing Mn concentration as a function of the depth 
could lead to an increase of D a £1 

Finally, we note, since B^q = in the simulations, 
the magnetization precesses circularly for ip = 90° and 
thus Im(m*m2 — mimj) = 2 sin 2 with the precession 
cone angle r. For all other orientations, m precesses 
elliptically which is accounted for in our simulations. In 
the simulations of the precession amplitudes, we have 
assumed an externally applied microwave magnetic field 
with [ioh = 0.1 mT. 

At an external field orientation of ip = 50° the uni- 
form resonance field is nearly constant across the layer, 
and consequently only one SWR mode is observed with 
an almost uniform magnetization precession across the 
layer, cf. Fig. [6](b). The precession amplitude is a mea- 
sure for the SWR intensity. While the fundamental mode 
at ip = 90° exhibits a larger precession cone at the in- 
terface, it rapidly decays as a function of the depth, in 
contrast to the nearly uniform precession amplitude for 
ip = 50°. Since the entire layer contributes to the power 
absorption, consequently, the SWR mode at ip = 50° is 
more intense than the fundamental mode for ip = 90°, 
which is indeed observed in the experiment [cf. Fig. [6] (d) 
and (e)]. 

For the magnetic field within the film plane [tp = 0°, 
cf. Fig.[6](c)], the uniform resonance field again varies lin- 
early across the film, however in a less pronounced way 
than for the out-of-plane field orientation and with an 
opposite sign of the slope. The sign reversal of the slope 
can be understood in terms of the uniaxial out-of-plane 
anisotropy parameter -Booi : positive values of these pa- 
rameters lead to an increase (decrease) of the resonance 
field for the magnetization oriented perpendicular (par- 
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TABLE I: Simulation parameters and their z dependence of the samples under study as obtained by fitting the simulations to 
the SWR measurements. For the anisotropy parameters the capital letters denote the value at the surface of the film and the 
lower case letters the slope as described in the text. For sample A, the first value of &001 was used for the first 100 nm and the 
second one for the remaining layer. In addition to the anisotropy parameters, the saturation magnetization is also assumed to 
vary linearly across the layer, while its absolute value is unknown and not important for the SWR simulations. 
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A 
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0.05 


25 


-0.3 


35 


0.09 


-3 


B 


1.3 


130 


-0.5 


-50 











20 


0.06 


-4 


C 


1.5 


75 


-0.4 


-55 


-0.04 


-15 





40 


0.11 


-4 


D 


3.5 


91 


-0.3 


-55 


-0.04 


-15 





20 


0.09 


-3 



allel) to the film plane, accounting for the slopes of the 
uniform resonance fields in Fig. [5] Since the gradient in 
the uniform resonance field is less pronounced for tp = 0° 
than for ip — 90°, the spin- wave modes are not resolved 
for ip = 0°, since their spacing is smaller than the SWR 
linewidth, leading to one rather broad line [cf. Fig. [B] (c) 
and (f)]. A steeper gradient in B411 in combination with 
a different Gilbert damping (or with an additional inho- 
mogeneous damping parameter) and amplitude scaling 
factor, could improve the agreement of simulation and 
experiment in the in-plane configuration, as discussed 
later. A detailed study of the in-plane anisotropy pro- 
file is however beyond the scope of this work. Given that 
the presented simulations were obtained with one set of 
parameters, the agreement of theory and experiment is 
reasonably good also for the in-plane configuration, since 
salient features of the SWR lineshape are reproduced in 
the simulation. 

Having discussed the angle- dependence of the SWR 
spectra, we turn to the z dependence of the out-of-plane 
anisotropy of sample A. Our simulations reveal that it 
is governed by the z dependence of both Bqqi(z) and 
Bi±{z). Assuming only a gradient in -B001 results in a 
reasonable agreement of theory and experiment for the 
external field oriented along [001] and [110], but fails to 
reproduce the spectra observed for the intermediate field 
orientations, e.g., ip = 50°. This is illustrated by the 
dashed black line in Fig. [5] (e), which represent simu- 
lations with a constant B^±(z) for ip = 50°. As can 
be seen, this simulation produces several spin-wave res- 
onances, whereas in the experiment only one resonance 
is present, which is better reproduced by the simulation 
with both Bqqi(z) and B4±(z) varying across the layer. 

We will now discuss the anisotropy parameters of all 
samples. In contrast to sample A, the out-of-plane 
anisotropy profile of all other samples appears to be gov- 
erned by a gradient in Bqqi(z). As already discussed 
qualitatively, the hard axis of the samples is determined 
by an interplay of -B001 an d B^±. For sample A and B 
B^± is positive and zero, leading to an out-of-plane hard 
axis. In contrast, sample C and D exhibit a negative B^±, 
leading to a hard axis between out-of-plane and in-plane. 
The £?4|| parameter is negative and of similar magnitude 
for all samples. 



Since the out-of-plane anisotropy profile of sample A 
is governed by Bqqi(z) and £?4j_(z), a comparison of the 
out-of-plane anisotropy profile between all samples based 
on anisotropy parameters is difficult. We therefore com- 
pare the uniform resonance fields, where both anisotropy 
parameters enter. As evident from Fig. |B1 the strongest 
influence of the magnetic inhomogeneity of the layers on 
the uniform resonance fields is observed for the exter- 
nal field along [001]. To compare the hole concentration 
profile in Fig. @] (a) with the anisotropy profile, we there- 
fore plot in Fig. [4] (b) the z dependence of the uniform 
resonance field hqH®®1 for this field orientation. The 
figure demonstrates that the gradient in /io-H^rri is cor- 
related with the gradient in p. For the sample with the 
strongest gradient in p the gradient in /j,qH®®1 is also 
most distinct while the samples with a weaker gradient 
in p exhibit a less pronounced gradient in fioH^ . How- 
ever, for sample D, exhibiting a nearly constant p, we 
still observe standing spin wave resonances for /j,qH\ | [001] 
[Fig. [11(d)], reflected in a slight gradient of fj, H°®l. This 
observation suggests that aditionally other mechanisms 
lead to a variation of the anisotropy profile. One possi- 
bility would be a gradient in the elastic strain of the layer, 
due to a non-homogeneous incorporation of Mn atoms in 
the lattice. However, x-ray diffraction measurements of 
this sample, in combination with a numerical simulation 
based on dynamic scattering theory, reveal a variation of 
the vertical strain Ae zz as small as 3 x 10~ 5 across the 
layer. According to the measurements in Ref. [T3 . such a 
variation in strain would lead to a variation of the -B001 
parameter by a few mT only, insufficient to account for 
the variation of /xo-Huni by almost 100 mT across the 
layer. A more likely explanation seems to be a varia- 
tion of the saturation magnetization, which should also 
influence the anisotropy parameters. In the simulation, 
a non-homogeneous saturation was assumed, potentially 
explaining also the observed gradient in the anisotropy 
parameters and therefore in the uniform resonance field. 

In contrast to the out-of-plane anisotropy parameters, 
S411 was found to depend only weakly on z, for all sam- 
ples except sample B where it was constant. Addition- 
ally, -B110, typically of the order of a few mT, 13 might 
have an influence and interplay with S411 in determining 
the in-plane anisotropy. We here however focus on the 
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FIG. 6: Simulated magnetization mode profile and uniform resonance field of sample A. The contour plots show the magneti- 
zation precession amplitude Im(mjm2 — mimj) as a function of the position within the film and the external magnetic field 
for the external field aligned (a) along [001], (b) at an angle of 50° with respect to [110] (cf. the inset in Fig. [5]) and (c) along 
[110]. The blue, dashed lines in (a)-(c) show the uniform resonance field, obtained by numerically solving Eq. © for each given 
field orientation. The dotted black lines in (a) indicate the resonance magnetic fields. In (d)-(f), a magnification of the data 
(blue dotted lines) and simulation (red solid lines) from Fig. 0(a) is shown using the same scale for all orientations. In (e), a 
simulation with a different set of parameters is shown for comparison (black, dashed line), see text. 
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out-of-plane anisotropy and therefore neglect -Biio in our 
simulations. An in-plane rotation of the external field 
would be required for a more accurate measurement of 
S 4 || and i?iio, but is outside the scope of this work. 

According to the valence-band model in Ref. [H, an 
oscillatory behavior of the magnetic anisotropy parame- 
ters is expected as a function of p. Therefore, depend- 
ing on the absolute value of p, different values for, e.g., 
dBoai/dp are expected. In particular, there are regions 
where a anisotropy parameter might be nearly indepen- 
dent of p and other regions with a very steep p depen- 
dence. Since the absolute value of p is unknown, a quan- 
titative discussion of the p dependence of the obtained 
anisotropy parameters based on the model in Ref. [12 
is not possible. In addition to p, the p-d exchange 
integral, 12 which may also vary as a function of the depth 
in a non-homogneous film, also influences the anisotropy 
parameters? 1 ^ further complicating a quantitative analy- 
sis. 

For all samples, we used a constant exchange stiffness 
D s in our modeling. As alluded to above, there is some 
ambiguity in this assumption, since the exchange stiff- 
ness and the gradient in the anisotropy both influence 
the mode spacing. For simplicity, however, we intended 
to keep as many simulation parameters as possible con- 
stant. The absolute values obtained for the exchange 
stiffness agree within a factor of 2 with the ones obtained 
in previous experiment s 36 ' 50 but are a factor of 2-4 larger 
than theoretically predicted.— For the reasons discussed 
above, there is a large uncertainty also in the derivation of 
the absolute value of D s from standing spin-wave modes 
in layers with a gradient in the magnetic anisotropy con- 
stants. 

In order to use one parameter set for all field- 
orientations, the Gilbert damping parameter was as- 
sumed to be isotropic in the simulations. The modeling 
of the SWR data could be further improved by assum- 
ing a non-isotropic damping, its value being larger for 
/io-H"||[110] than for ^ H\\[001} [cf. Fig. 0. This how- 
ever, only improves the result when assuming a field 
orientation-dependent scaling factor for the amplitude, 
which could be motivated, e.g., by the assumption that 
the microwave magnetic field present at the sample po- 
sition depends on the sample orientation within the 
cavity. The absolute values of a determined here are 
comparable with the ones obtained by ultra-fast opti- 
cal experiments* 5 - 2 - but are larger than the typical a — 
0.01. ..0.03 values found by frequency-dependent FMR 
studies . 53 ' 54 As already alluded to, inhomogeneous line- 
broadening mechanisms may play a dominant role^ in 
particular for as-grown samplesi 5 - 5 - We therefore assume 
that the values for a obtained in this study overesti- 
mate the actual intrinsic Gilbert damping. A frequency- 
dependent SWR study would be required to determine 
the intrinsic a. Such a study could possibly also reveal a 
p-dependent a as theoretically predicted.— In our study, 
assuming a z dependent a did not improve the agree- 
ment between simulation and experiment, corroborating 



the conjecture that inhomogeneous broadening mecha- 
nisms dominate the linewidth and therefore obscure a 
possible z dependence of a. 



IV. SUMMARY 

We have presented a finite difference- type modeling ap- 
proach for standing spin- wave resonances based on a nu- 
merical solution of the LLG equation. With this generic 
formalism, SWR spectra can be simulated accounting for 
elliptical magnetization precession, for arbitrary orienta- 
tions of the external magnetic field, and for arbitrary pro- 
files of all magnetic properties, including anisotropy pa- 
rameters, exchange stiffness, Gilbert damping, and sat- 
uration magnetization. The approach is applicable not 
only to (Ga,Mn)As but to all ferromagnets. 

Four (Ga,Mn)As samples, epitaxially grown with V/III 
flux ratios of 1.1, 1.3, 1.5, and 3.5 were investigated by 
ECV and spin-wave resonance spectroscopy, revealing a 
correlation of a linear gradient in the hole concentration 
with the occurrence of standing spin wave resonances, in 
particular for the external field oriented out-of-plane. Us- 
ing the presented modeling approach, the SWR spectra 
could be reproduced in a simulation with one parameter 
set for all external field orientations. The simulation re- 
sults demonstrate that the profile of the out-of-plane uni- 
form resonance field is correlated with the hole concentra- 
tion profile. However, our measurements and simulations 
show, that a non-uniform hole concentration profile is not 
the only cause that leads to the observed non-uniform 
magnetic anisotropy; possibly, a variation in the satura- 
tion magnetization also influences the anisotropy param- 
eters. To gain a quantitative understanding of this issue, 
more samples with known hole concentrations would be 
required, where both the absolute values and the profiles 
of p are varied. Such a study was, however, outside the 
scope of this work. 

Besides the modeling of SWR intensities and 
linewidths, the presented formalism yields the magne- 
tization precession amplitude as a function of the po- 
sition within the ferromagnet. It can therefore be used 
to investigate spin-pumping intensities in (Ga,Mn)As/Pt 
bilayers.— The spin-pumping signal, detected as a volt- 
age across the Pt layer, should be proportional to 
the magnetization precession cone in the vicinity of 
the (Ga,Mn)As/Pt interface. By measuring the spin- 
pumping signal as well as the SWR intensities of 
(Ga,Mn)As/Pt and by using our modeling approach, it 
should be possible to investigate to which extent a mag- 
netization mode which is localized at a certain posi- 
tion within the (Ga,Mn)As layer contributes to the spin- 
pumping signal. 
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Appendix A: Coordinate Transformation and Free 
Enthalpy derivatives 

The transformation between the crystallographic coor- 
dinate system [x, y, z) and the equilibrium system (1,2,3) 
is given by 



m x \ 




(mi \ 




m y 


= T 


m 2 , 


(Al) 


m z 1 




[m 3 l 





with 



T 



COS t) COS 00 

cos 9q sin 0o 
— sin 6>o 



— sin 0o sin uq cos 0o 
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(A2) 



The derivatives of the free enthalpy density Eq. © with 

respect to the magnetization components are 

G3 — 9 m3 G\ m=mo = — HqH 3 + 2B 00 i cos 2 6 

+ B x jo (sin 9o cos 0o — sin #0 sin 0o ) 2 

+ 4B 4± cos 4 8 

+ 4B 4 [| sin 4 6» (cos 4 O + sin 4 O ) (A3) 

G2I = G12 = d mi dm 2 G\m=-m 

— cos ^o(l ~ 2 cos 2 0o) [-Bjjo 

+ 12B 4 || sin 2 6»o cos O sin O ] (A4) 
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+ Bi 10 cos 2 6*o (cos 0^o — sin 0o ) 2 

G22 = d m2 d m2 G\ m=mo = 2_B llQ (sin0^o + cos0 o ) 

+ 24B411 sin 2 6»o cos 2 O sin 2 ■ (A6) 
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In this Appendix, we describe how the spin- wave equa- 
tion can be numerically solved by the finite difference 
method. We start with the simple case of a circulary pre- 
cessing magnetization, neglecting Gilbert damping and 
the driving field (Sec. IB ip . Then we turn to the gen- 
eral case, where the magnetization precesses elliptically 
and the Gilbert damping as well as the driving field are 
included fSec. lB2j) . 

1. The One-Dimensional, Homogeneous, 
Undamped Case 



Here, we describe how the resonance fields and the 
spin-wave modes can be found, assuming a circularly 
precessing magnetization m,2 = im\ = m, a constant ex- 
change stiffness, and a z independent equilibrium magne- 
tization. This case has been considered in Ref. |36| using a 
semi-analytical approach to solve the spin- wave equation 
Eq. (|10p . The approach considered here, is slightly more 
general, as it is straight forward to determine resonance 
fields and eigenmodes of the system for an arbitrary z 
dependence of the uniform resonance field. To solve 
Eq. (TiT)]) , we divide the ferromagnetic film into a finite 
number n of layers with equal thickness I and constant 
magnetic properties within each of these layers. The z 
dependence of m and /JoH^] is thus given by an index 
j = 1...7J. Within each of these layers the uniform reso- 
nance field and m(z) are thus constant and given by the 
values ^oH uni 3 =: iP and m J , respectively. The second 
derivative of rh is approximated by 



m" (z = j ■ I) 



p- 1 - 2m? + m J+1 



(Bl) 



Consequently, Eq. (flOf is converted to the homogeneous 
equation system 



KJ- 1 + 2d -d 
-d K j +2d -d 
-d Ki +1 + 2d 



( : \ 

m?~ 
m? 
m- J+1 



fJ-o 



H 



m? 



(B2) 



with the abbreviation d = —D s /l 2 . The boundary condi- dition) reads as m° = m 1 and rh n 1 = m n and can be 
tion of natural freedom^ (von Neumann boundary con- incorporated in Eq. (IB2[) . Since the matrix on the left 
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hand side of Eq. (|B2[) is sparse, it can be efficiently diag- 
onalized numerically, yielding the resonance fields (eigen- 
values) and the corresponding modes (eigenvectors). Af- 
ter diagonalizing the matrix, the relevant resonance fields 
are found by sorting the eigenvalues and considering only 
the modes with positive resonance fields, corresponding 
to the bound states in the particle-in-a-box analogon. 
The SWR amplitude of each mode is proportional to its 
net magnetic moment; thus, the amplitudes can be found 
by integrating the (normalized) eigenmodes. The mode 
profile, the resonance fields, and the SWR intensities are 
illustrated in Fig. [5] for a constant and a linearly varying 
uniform resonance field. The finite linewidth of the SWR 
modes can be accounted for by assuming a Lorentzian 
lineshape for each mode with a certain linewidth and 
with the resonance fields and intensities calculated as 
described above 35 . Note that this approach to derive 
resonance fields and intensities is only valid if the mode 
separation is large compared with the linewidth of the 
modes; this restriction does not apply to the model pre- 
sented in the Appendix IB 21 



netic film into a finite number n of layers with equal 
thickness I and constant magnetic properties within each 
of these layers. In contrast to the case in the Appendix 
IB 11 where only the uniform resonance field was varied 
across the layer, here potentially all magnetic proper- 
ties entering Eq. J7| can be assumed to be z dependent. 
Additionally, the components of the driving field fJ,ohi 
(i = 1,2), can also vary as a function of z, since the 
(1,2,3) frame of reference is z dependent and thus the 
projections of the driving field have to be calculated for 
each layer. The z dependence of the components mi 
(i = 1,2), of the parameters H\\, ffia, H21, #22 (de- 
fined in Sec. [TTJ) and the exchange stiffness is thus given 
by the index j = 0...n; the second derivative of each of 
the components nrii is approximated as in Eq. (|B1|) . 



2. The General Case 

To solve Eq. ([7]) for arbitrary /j.qH and arbitrarily vary- The linearized LLG equation Eq. ([7]), is thus converted 

ing magnetic properties, we again divide the ferromag- into the inhomogeneous equation system 
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with the abbreviation d? = —D(/l 2 . At the boundaries 
of the magnetic film we again assume the spins to exhibit 
natural freedom mj = m\ and m" = 

To simulate a spin-wave spectrum for a given orienta- 
tion of the external field and a given profile of the mag- 
netic properties, we numerically sweep the magnetic field 



and calculate the equilibrium magnetization orientation 
for all indices j — 0...n at a given external field. The 
inverse of the matrix in Eq. (|F33|) . multiplied by hqM(z), 
is the generalized Polder susceptibility tensor xit^oH , z), 
which relates the transverse magnetization with the driv- 
ing field, cf. Eq. (fTTj) . 
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